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ABSTRACT 

The stellar kinematics of the spheroids and discs of SO galaxies contain clues to 
p^ . their formation histories. Unfortunately, it is difficult to disentangle the two compo- 

QQ ■ nents and to recover their stellar kinematics in the faint outer parts of the galaxies 

fT) \ using conventional absorption line spectroscopy. This paper therefore presents the 

stellar kinematics of six SO galaxies derived from observations of planetary nebulae 
(PNe), obtained using the Planetary Nebula Spectrograph. To separate the kinemat- 
^^ I ics of the two components, we use a maximum-likelihood method that combines the 

^^r) ■ discrete kinematic data with a photometric component decomposition. The results of 

this analysis reveal that: the discs of SO galaxies are rotationally supported; however, 
the amount of random motion in these discs is systematically higher than in compa- 
rable spiral galaxies; and the SOs lie around one magnitude below the Tully-Fisher 
^^ I relation for spiral galaxies, while their spheroids lie nearly one magnitude above the 

%^ . Faber- Jackson relation for ellipticals. All of these findings are consistent with a sce- 

_ C^_' nario in which spirals are converted into SOs through a process of mild harassment or 

"pestering," with their discs somewhat heated and their spheroid somewhat enhanced 
by the conversion process. In such a scenario, one might expect the properties of SOs 
to depend on environment. We do not see such an effect in this fairly small sample, 
although any differences would be diluted by the fact that the current location does 
not necessarily reflect the environment in which the transformation occurred. Similar 
observations of larger samples probing a broader range of environments, coupled with 
more detailed modelling of the transformation process to match the wide range of 
parameters that we have shown can now be measured, should take us from these first 
steps to the definitive answer as to how SO galaxies form. 

Key words: galaxies: elliptical and lenticular - galaxies: evolution - galaxies: kine- 
matics and dynamics. 
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1 INTRODUCTION 

The origin of lenticular, or SO, galaxies, and in particular 
whether they are more closely related to spiral galaxies or 
elliptical galaxies, remains obscure. They display the bulge- 
plus-disc morphology that we associate with spiral galaxies, 
but they lack the young stars in spiral arms that we asso- 
ciate with such systems, and tend to be more dominated by 
their spheroidal bulge components, making it tempting to 
associate them with elliptical galaxies. 

Moving away from purely morphological considerations, 
we might hope that kinematic information can help to ascer- 
tain which other galaxies are SOs closest cousins, since the 
traditional view is that disc-dominated spiral galaxies will 
reflect that morphology in rotationally-dominated kinemat- 
ics, while elliptical galaxies are largely supported by ran- 
dom motions. However, more recent data suggests that the 
situation is some what subtler. Th e impo rtance of rotation 
was quantified by lEmsellem et al\ (|2007|) , who defined the 
quantity \r^ the cumulative specific angular momentum of 
galaxies as a function of radius, seeking to draw a distinc- 
tion between disc- and bulge-dominated systems via this 
parameter. Unfortunately, they show that there is no clear 
dichotomy in this quantity, with a broad spectr um of de- 
grees of rotational support. ICoccato et al\ (| 20091 ) calculate 
the same quantity for a sample of elliptical galaxies, but us- 
ing PNe as tracers of the velocity field and, in this way, ex- 
ploring a wider range of radii. They show that the recovered 
\r smoothly join the quantity obtained with the SAURON 
data, which describe the central regions of the galaxies, 
where PNe are not detectable. This radially-extended recon- 
struction of the specific angul ar momen tum for a sample o f 
elliptical (|Coccato et a/.ll2009l ) and SO (|Cortesi et a/.ll201lh 
galaxies, shown in Figure [U displays the same continuous 
range of profiles, seemingly unrelated to the galaxy morphol- 
ogy. It also presents a surprising degree of variation in profile 
shape at the large radii probed by PNe, with some drop- 
ping dramatically, apparently changing from rotationally- 
supported systems to being dominated by random motions 
in their outermost parts. 

We first encountered this phenomenon in an anal- 
ysis of the SO galaxy NGC 1023, where a pilot study 
of PNe with the Planetary Nebula Spectrograph (PN.S) 
showed a major decrease in rotational motion at large radii 
((Noordermeer et al. 2008). This drop in angular momen- 
tum was initially interpreted as evidence that a normal 
rotationally-supported stellar disc had been dramatically 
heated into random motion by a merger in its outer parts. 
However, by more careful modelling of separate bulge and 
disc components, we were able to show that this behaviour 
was in fact caused simply by the superposition of a cold ro- 
tating disc and a hot spheroidal component, whose relative 
contributions to the observed kinemat ics varied with bot h 
radius and azimuth around the galaxy (jCortesi et a/.ll201lh . 

This discovery underlined the importance of treat- 
ing SO galaxies as multi-component systems (see also 
iKormendv fc Bender! (i2012)), and taking care to separate 
out the distinct stellar bulge and disc elements kinemati- 
cally as well as photometrically. Once this decomposition 
has been performed properly, we have access to a range of 
diagnostics that can be used to provide clues to the forma- 
tion history of these galaxies. For example, if they are sim- 
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Figure 1. Plot showing a measure of specific angular momentum 
as a function of radius for a sample of early- type galaxies. Black 
solid lines ar e SOs from | Cortesi et oU (|201lh , black dashed lines 
are SOs from ICoccato et al\ (|2009l V and red lines are ellipticals 
from lCoccato etoU (|2009l ). The dotted line shows the suggested 
separation bet ween slow and fast rot ators as a kinematic classifi- 
cation scheme (|Emsellem et aLll2007T ). 



ply spiral galaxies that have quietly ceased forming stars, 
one would expect their stellar discs to have the same kine- 
matic properties as those in spirals; if, on the other hand, a 
more violent event such as a minor merger led to the tran- 
sition from spiral to SO, one might expect the disc, if not 
destroyed in the process, to at least hav e had its random 
motions increased significantly (Bournau d, Jog k, CombesI 
I2OO5I ). Equally, with reliable kinematic parameters for the 
individual components, we can see where they lie relative to 
the usual kinematic scaling relations followed by spirals and 
ellipticals, to see if they could plausibly have evolved from 
such progenitors. 

In this paper, we therefore apply the kin ematic bulge- 
disc d ecomposition technique developed in ICortesi et al\ 
(|201l[ ) to the sample of six SO s with suitable PNe data pre- 
sented in Cortesi et al. (2013 j. This sample was selected to 
span a range of environments in order to see if the evolu- 
tionary properties depend systematically on surroundings: 
NGC 7457 and NGC 3115 are isolated; NGC 1023 and 
NGC 2768 are the dominant galaxies of two small groups 
(5-6 members); while NGC 3384 and NGC 3489 are satellite 
galaxies of the Leo Group (30 members). The Hubble types 
of these galaxies and the adopted distance moduli are listed 
in Table [T] 

The remainder of this paper is laid out as follows. In 
Section [21 we present the photometric bulge-disc decompo- 
sition that underlies the analysis, and confirms the nature of 
these galaxies as composite SO systems. Section [3] then uses 
this decomposition to model the spheroid and disc kinemat- 
ics of each system, and Section |4] looks at the characteristic 
kinematics of the separate components for indications as 
to their origins. Section [5] brings this evidence together to 
present the emerging picture as to the steps leading to the 
formation of a lenticular galaxy. 
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2 PHOTOMETRIC DECOMPOSITION 

In order to correctly assign each PN's observed velocity to 
disc or bulge kinematics, we must first calculate the rela- 
tive contributions of disc and bulge to the total light at its 
location. To do so, we follow t he same photometric decom- 
position procedure described in lCortesi et al\ (|201lh , briefly 
summarised here for completeness. 

Where possible, the galaxies' images used for the de- 
composition are K-h and data from the 2MASS survey 
([Skrutskie et a/.l2006l ) , since these near infrared images min- 
imize the influence of any modest amounts of dust extinc- 
tion present in these SO systems. For NGC 3489, the 2MASS 
image did not go deep enough to carry out the decomposi- 
tion, so we instead used the reddest SDSS optical image, 
and for NGC 1023 w e also employed a deep R-band image 
([Noordermeer et a/.l 12008). which allowed us to mask the 
compani on galaxy NGC 1023 A while performing the decom- 
position ([Cortesi et a/.ll201lh . The details of the images and 
other input data used in the fitting process are presented in 
Table m 

The images were modeled using GALFIT ( Peng et al\ 

^2002) to fit an exponential disc and a Sersic profile spheroid. 
For most galaxies, a satisfactory fit was obtained by fixing 
the Sersic profile index at the conventional de Vaucouleurs 
law value of n = 4, but in the case of NGC 2768 this pro- 
duced a poor fit, so the index was left as a free parameter. 
The resulting fits to the images are presented in Figure [21 
and the associated values for the best fit parameters are 
given in Table [2] As is apparent from Figure O in some 
cases this simple two component model fits the galaxy very 
well, while in others systematic residuals indicate that the 
system is somewhat more complex. In particular, we find: 



• NGC 3115 contains an additional very thin disc- like 
structure as well as the thicker disc component that we have 
fitted. In addition, the de Vaucouleurs fit is not perfect for 
the bulge at very small radii. 

• NGC 3489 contains a faint but significant ring struc- 
ture. 

• NGC 3384 shows a central dipolar structure in the resid- 
uals, seemingly indicative of an off-centre nucleus. 

Although all these features are interesting, and tell us that 
even the plainest looking SO galaxy can be quite complex, 
they are all either localized in regions of high surface bright- 
ness where we do not detect the PNe used in this kinematic 
analysis, or they are of low surface brightness compared to 
the main bulge and disc components. Accordingly, they do 
not compromise our ability to use the simple two compo- 
nent fit to determine the relative contributions of disc and 
spheroid to the light at each point in the galaxy. 

To illustrate the assignment of PNe to the two com- 
ponents. Figure [3] shows a greyscale image of the spheroid- 
to-total light at each point in the model. The PNe are also 
plotted and identified by whether they lie in a region where 
spheroid or disc dominates, but note that, in the kinematic 
likelihood analysis below, each is assigned an exact value 
fi from this image, such that it has a probability fi that 
it comes from the spheroidal co mponent and 1 — fj that it 
comes from the disc component ([Cortesi et a/.ll201lh . 



3 DISC AND SPHEROID KINEMATIC 
PROFILES 

Having obtained the disc-spheroid light decomposition for 
each galaxy, we can now fit all the line- of- sight velocities 
for the PNe shown in Figure [3] using a maximum likelihoo d 
fit. The method, described in detail in lCortesi et al\ (|201lh , 
essentially involves fitting in radial bins using a model com- 
prising a simple Gaussian spheroid line-of-sight velocity dis- 
tribution plus a similar Gaussian velocity distribution for 
the disc component. For the disc component, we allow a 
non-zero mean velocity to fit rotation, varying with azimuth 
as geometrically required for a rotating disc, and also incor- 
porate the fact that its line-of-sight velocity dispersion is a 
different projection of the disc's coupled radial and tangen- 
tial velocity dispersions, ar and ad) respectively, at different 
azimuths. We neglect the contribution of the z-component 
of the velocity dispersion: it is intrinsically smaller, and its 
modest projection along the line of sight further reduces its 
significance in these inclined systems. We also assume that 
the disc is not too hot, so we can invoke the epicycle ap- 
proximation to couple the values of ar and ad). This fitting 
process thus solves for a simplified model of both disc and 
spheroid kinematics, and also allows us to identify and re- 
ject PNe that do not fit the model, based on their individual 
contributions to the total likelihood. 

Any PNe rejected in this process are highlighted 
in Figure [3l and the final resulting best-fit kinematic 
parameters as a function of radius are shown in Fig- 
ure m In each case, the radial bins have been chosen to 
ensure that each contains at least 30 PN e , as f ound to 
be a suitable minimum in ICortesi et all (|201lh , which 
is why the bin sizes and number vary from galaxy to 
galaxy. The resulting elliptical bin boundaries are shown 
in Figure (3] For comparison. Figure [H also shows conven- 
tional absorption- line kinematics along the major axes 
from various sourc e s (Caon. Macchetto fc Pastoriza 200d: 
Simien fc Prugniell 1 19971: iNorris, Sharpies fc KuntschneJ 



20061 : toebattis ta. Corsini fc Aguerril l2002l ). Since the disc 
light dominates in this region, we compare these data to the 
derived disc kinema tics. In one case there is also published 
minor-axis data ([Norris. Sharpies fc Kuntschneil |2006|) 
where the bulge light dominates, so the absorption-line 
dispersion profile is compared to the bulge kinematics that 
we derive. In all cases, the agreement between the two 
methods is good, but the comparison underlines how con- 
ventional absorption-line spectroscopy is typically limited 
to the bright inner parts of each galaxy. 

One simplification that we have made in this analy- 
sis is in assuming that the spheroidal component is non- 
rotating. In general, there is not sufficient data to allow 
us to relax t his assumption, but in the case of NGC 2768 
iForbes et al\ ([2012) found that the larger number of PNe 
allowed this extra degree of freedom to be introduced. How- 
ever, the resulting spheroidal component was found to be 
completely dominated by random motions, so the differ- 
ence was negligible, justifying the assumption in the case 
of this galaxy and rendering it plausible for the other galax- 
ies where there were not sufficient data to test it directly. 
This assumption also fits with the findings in spiral galax- 
ies, where classical bulges tend to be very slowly rotat- 
ing 



where classical bulges tend to be very slowly rotat- 
([MacArthur, Gonzalez fc Courteaull2009l l while pseudo- 
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Name 


Type 


Distance modulus 

[mag] 


Archive 


Band 


Angular scale 
[" per pixel] 


Zero point 

[mag] 


NGC 3115 


SO-edge-on 


29.77 


2MASS 


K 


1 


20.86 


NGC 7457 


SA(rs)0- 


30.45 


2MASS 


K 


1 


20.09 


NGC 2768 


E6 


31.59 


2MASS 


K 


1 


20.10 


NGC 1023 


SB(rs)0- 


30.13 


2MASS 


K 


1 


20.03 


NGC 3489 


SAB(rs)0+ 


30.25 


SDSS 


z 


0.4 


22.55 


NGC 3384 


SB(s)0- 


30.16 


2MASS 


K 


1 


20.09 



Table 1. Basic data on the sample and images analysed. From left to right: galaxy name, galaxy type, distance modulus from lTonrv et al\ 
(|200lh shifted by 0.16 magnitudes (see text for details), archive, band, angular scale, zero point 
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NGC 3115 


8.34 


53.69 


0.39 


67 


45.00 


7.17 


26.19 


4 


0.31 


45.00 


0.74 


NGC 7457 


8.56 


27.07 


0.48 


62 


-57.28 


9.49 


11.62 


4 


0.62 


-46.04 


0.30 


NGC 2768 


8.19 


42.93 


0.29 


73 


-86.25 


7.23 


50.46 


4.65 


0.66 


-85.39 


0.71 


NGC 1023 


7.02 


59.08 


0.26 


74 


84.12 


6.9 


17.86 


4 


0.75 


75.59 


0.53 


NGC 3489 


8.18 


24.98 


0.49 


61 


-15.05 


8.05 


7.64 


4 


0.72 


-19.50 


0.54 


NGC 3384 


7.92 


63.73 


0.34 


70 


52.50 


7.29 


15.20 


4 


0.83 


60.51 


0.64 



Table 2. Results from GALFIT fit. From left to right: galaxy name [1], disc apparent magnitude [2], disc scale length [3], disc axes ratio 
[4], used to obtain the galaxy inclination [5], disc position angle [6], spheroid apparent magnitude [7], effective radius [8], Sersic index, 
equal to 4 where this value fitted well, [9], spheroid axes ratio [10], spheroid position angle [11] and bulge-to-total light ratio [12]. For 
NGC 3489 the recovered magnitudes in z-band have been colour corrected. 



bulge s with low Sersic ind ices display more rotational sup- 
port (|Fabricius et a/.ll2012l ) ; the high Sersic indices for these 
galaxies (see Table [2} suggests we might expect a corre- 
sponding lack of rotation. 

One further assumption we made is that the galaxies 
have only two components, a disc and a spheroid. As dis- 
cussed above, the images show some departures from this 
simple model, although in most cases these extra features 
are fairly minor. The only exception is NGC 3115, which 
seems to contain two discs of different thickness. This object 
would clearly be interesting to study in its own right with 
a more sophisticated model, but for the current analysis, it 
makes more sense to treat it in a manner consistent with 
the other galaxies, fitting only the spheroid and the brighter 
disc. It is interesting to note that this model results in more 
PNe being rejected from the fit than in other galaxies, which 
once again illustrates the power of the likelihood approach 
to reject discrepant objects and home in on the c omponents 
being fitted, as was originally seen in lCortesi et a l (2011). 

As is apparent from Figure [H this analysis reveals a re- 
markably consistent kinematic picture amongst these galax- 
ies. The one exception seems to be NGC 3489, although 
it is notable that this galaxy contains the smallest number 
of available kinematic tracers, which may be compromis- 
ing the results somewhat. However, apart from this galaxy, 
all the systems here display kinematics with many features 
in common: the spheroidal components have a dispersion 
that varies very little with radius, and the disc components 
rise to flat rotation curves with ordered motions dominat- 
ing over random velocities at all radii. It is interesting to 
compare this situation with the much more heterogeneous 
picture presented in Figure [Tl it would appear that at least 
some of the variety in that plot arose from the net effect 
of superimposing multiple kinematic components whose re- 
spective contributions vary from galaxy to galaxy. This new 
consistency provides some confidence in the analysis, but 



also offers at least the possibility that these galaxies are 
sufficiently homogeneous for some common underlying for- 
mation mechanism to exist. 



4 ANALYZING THE CHARACTERISTIC 
KINEMATICS 

4.1 Deriving characteristic values 

Having determined the kinematic profiles of the spheroid 
and disc components of these SO galaxies, we can now start 
to use these data to seek archaeological evidence as to how 
they formed. As a starting point, we translate these profiles 
into characteristic values for the kinematics of each compo- 
nent. For the disc, we determine the mean streaming mo- 
tion, K, and the two components of velocity dispersion, ad) 
and ar, at three disc scale lengths as determined by the 
photometric parameters (see Table [2]) , since by this radius 
they seem to have settled to their asymptotic values. The 
only exception is NGC 3489, where the limited amount of 
kinematic data restricts us to calculating these quantities 
at 2Rd; in practice this makes little difference, as in the 
other galaxies the parameters do not change significantly 
over this radial range. The resulting values are presented in 
Table O In the spheroids, the velocity dispersion does not 
change significantly with radius, so we simply calculate a 
luminosity- weighted average value. 



^2 _ / jj ^sph^CTsph.i 
^ sph 



E.i 



sph,i 



(1) 



where Lsph,i is the luminosity of the spheroid in each radial 
bin, as ascertained from the photometric fit. The resulting 
values of dsph are also listed in Table (3] 

The other kinematic quantity we need to derive is the 
characteristic circular speed of each galaxy, 14, which pro- 
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Figure 2. GALFIT analysis of the SO sample. Left panels: galaxy image; middle panels: GALFIT model image; right panels: difference 
image, scaled to emphasize any residuals in the fit. The ratio between the flux in the residual image and in the scientific image is around 
5%. From top to bottom: NGC 3115 in a [800'' x 800''] box, NGC 7457 in a [600" x 600"] box, NGC 2768 in a [800" x 800"] box, 
NGC 3489 in a [~ 300" x 300"] box and NGC 3384 in a [c^ 600" x 600"] box. The dimension of the box corresponds to the dimension 
of the image in which the fit was performed, apart from NGC 3384, which was fitted simultaneously with the two companion galaxies 
appearing in the 2MASS [1200" x 1200"] image, in order to remove any residual contamination. 



vides a measure of the system's mass. Because of the pres- 
ence of significant random motions, we cannot simply use 
the mean streaming speed of the stars, but must correct 
these motions for asymmetric drift via the equation 



V^ 



K' + ctI 



-crr(l+ -71 + -71 ) 

a in r a in r 



(2) 



(jBinnev fc TremaindllQSTI ). If the random motions are not 
too large and the rotation curve is close to flat, then the 
epicycle approximation implies that ar — V2(j(j). If we fur- 
ther use the photometrically-derived exponential disc pro- 
file (see Section [2]) to determine z/(r) and fit the observed 
variation in velocity dispersion with radius using a further 
exponential with its own scalelength, 
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Figure 3. Probability map that a PN belongs to the disc (darker areas) or spheroid (lighter areas). The locations of detected PNe are 
also shown, with those more likely to belong to the disc shown as triangles, and those more likely to belong to the spheroid as squares. 
Colours indicate the line-of-sight velocities of the PNe. Open circles are drawn around PNe rejected by the likelihood fit to the kinematics. 
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Figure 4. disc and spheroid kinematics of the sample galaxies assuming a spheroid+disc model. For each galaxy, we show the rotation 
velocity of the disc (top panel), its coupled tangential, cr^^, and radial, dr dispersions (middle panel), and the velocity dispersion 
of the spheroid (bottom panel). Vertical error bars indicate uncertainty, while horizontal error bars show the radial b inning. Where 
available, we also s how kinematics derived from conventional absorption-line spectra (open circles) ( Caon, Macchetto Sz Pastoriza 201 



ISimien fc Prugniel|[l997l : iNorris, Sharpies fc KuntschnedlJOOel : iDebattista, Corsini fc AguerrillJOOJ 
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Name 


(K)disc 


iy^/(yct>) disc 


(Vc)disc 


f^sph 




[km/s] 




[km/s] 


[km/s] 


NGC 3115 


220+^9 


3.3±^o-^, 


264^g{ 


1071^^ 


NGC 7457 


113+^^ 


q 7+0.6 


136+^g 


54l«^ 


NGC 2768 


232+39 


q q + 0.6 


q-|f^+79 


123+11 


NGC 1023 


2AAt% 


^ q + 0.8 


274+^^ 


112lt^ 


NGC 3489 


144+g 


4.2+g-? 


171+^8 


571?? 


NGC 3384 


179+^28 


^ q+0.6 


196+^7 


89-19 



Table 3. Kinematic results. From left to right: galaxy name [1], rotation velocity in the disc at ?fRd [2], ratio between rotation and 
velocity dispersion along the tangential direction for the galaxy disc calculated at 2>R(i [3], circular speed in the disc at 37?^ [4], light 
weighted spheroid dispersion velocity [4]. 
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Figure 5. Ratio of stellar rotation velocity in the disc to the 
disc velocity dispersion for the SO galaxies of this sample (filled), 
and the corresponding quantity for a comparison sample of spiral 
galaxies from the literature (diagonally shaded) (Bo ttemal 119931 : 
iHerrmann fc Ciaxdullo 2009i), normalised to the same number as 
the SO galaxies. 



al = cr?(0)exp(---), 

the asymmetric drift equation simplifies to 



V^{r) 






(3) 



(4) 



We can hence estimate Vc at any radius using this equa- 
tion; Table [3] lists the derived characteristic value of this 
quantity for each galaxy at 3i?d, which matches the fiducial 
radius used for the other kinematic parameters; it is also 
at large enough radii that the circular speed will have con- 
verged to its characteristic asymptotic value. A check on the 
validity of this simplified equa tion i s provi ded by comparing 
our results to those of Davis et al\ (1201 ih, who carried out 



et g/jj 
ling ([C 



full anisotropic Jeans modelling (|CapDellaril l2008l ) on two 
of the sample galaxies. They obtained circular velocities of 
310 km s-^ for NGC 2768 and 160 km s"^ for NGC 3489, in 
good agreement with the values derived here. 



4.2 Disc kinematics 

As previously mentioned, one key diagnostic of the evo- 
lutionary past of a disc is provided by the ratio between 
its ordered and random motions: if they are simply spiral 
galaxies that have ceased forming stars, then they should 
be dominated by mean streaming motions in the same way 
as their progenitors, whereas any more violent transition 
such as one precipitated by a merger would tend to heat 
the disc and hence decrease the dominance of the rotation. 
Table [3] lists the characteristic values of V^/acf) for these SO 
systems, yielding a mean value of 4.2, with an RMS scatter 
of 0.8. These values con trast with the V^/acf) ~ 1 obtained 
bv lBournaud, Jog fc Co mbes (2005|) in the simulations that 
formed SOs from a range of minor mergers, suggesting that 
this is not the mechanism responsible for the creation of 
these galaxies. 

To determine whether the data are consistent with the 
alternative hypothesis that these discs are simply "dead" 
spiral discs, we consi der the k inema tics of eight spiral 
galaxies prese nted by JBottemal ([1993 ) and five shown by 
Herrmann & C iardulld (|2009r ), which together form a com- 
parison sample with a range of rotation velocities similar to 
the SOs presented here. These comparison data a re unfor- 
tunate ly somewhat heterogeneous in nature. The iBottemal 
(| 19931 ) kinematic values were derived at one disc scale length; 
since the velocity dispersion tends to fall with radius while 
rotation increases, this smaller fiducial radius will tend to 
produce smaller values of V^/acf) than at our adopted ra- 
dius of 3i^(j- The comparison is mo re direct with the data 
from Herrmann & Ciardulld (|2009|), as we can interpolate 
a value at ?>Rd from the kinematic profiles they present. In 
both cases, we again use the epicycle approximation to in- 
terchange between Gr and Gcf). Figure [5] shows the resulting 
estimates of V^/gcj) for both the spiral comparison sample 
and the current sample of SOs. Although there is a reason- 
able degree of overlap, implying that these SOs could have 
formed passively from spirals, there are also strong indica- 
tions that the values are systematically lower for the SOs, 
which is not the sense one would expect if it arose from the 
possible bias in the smaller fiducial radii used by Bott emal 
(1993). Thus, if these SOs formed from a random selection 
of comparable present-day spirals, it would seem that some 
heating of their discs must have occurred during the transi- 
tion. 

Although the sample becomes small when divided in 
this way, there is no evidence that the degree of heating 
of the disc depends on environments, since the SOs from 
isolated, small group and large group surroundings (as de- 
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Figure 6. X-band Tully-Fisher relation for the SO galaxies an- 
alyzed here, compared to the relation for spiral galaxies (blue 
circles and line of b est fit) from ( Tullv fc Piercj [2000|) and 
iRothberg et al\ (|2000[ ) 



scribed in Section [T]) do not display systematically different 
values of V^Igc^. We are thus left with a developing pic- 
ture in which the formation of an SO from a spiral appears 
somewhat more violent than a simple shutdown of star for- 
mation, but not as extreme as a merger-induced transition, 
with no indication that the mechanism depends strongly on 
environment. 



4.3 Scaling relations 

We can seek further insight into the possible scenario for 
transformation, and also try to address the question of 
whether SOs are more closely related to spirals or ellipti- 
cals, by looking at where these systems lie on the scaling 
relations respected by the other galaxy classes. 



Ji..S.l Tully-Fisher relation 

The simplest test we can carry out is to place these SO 
galaxies on the Tully-Fisher relation, which displays a very 
tight correlation between a spiral galaxy's lum inosity and 
its circular rotation speed (JTullv fc Fisherll 19771 ). Such com- 
parisons have previously been carried out by v arious au- 
thors, with somewhat mixed results (I Ar agon- Salamanca 
20081 : JNeistein et~al\ Il999l: IWilhams. Bureau &: Cappellaril 



;ne reason 



201ol : fpahre. Diorgovski fc de Carvalhol Il998l ). Or 

for the discrepant results is the difficulty in obtaining the 

circular rotation speed: for the spiral galaxies, this quantity 



can be reliably estimated from gas kinematics, but for the 
SOs we are reliant on stellar motions. The previous studies 
using conventional absorption-line spectroscopy have there- 
fore faced two significant problems. First, unlike PN data, 
the absorption-line spectra do not typically reach particu- 
larly large radii, so probe regions where the rotation curve 
has not reached its asymptotic value and where the random 
motions are still large so the somewhat-uncertain asymmet- 
ric drift correction is also large (see Section [4. 1|) . Second, as 
we saw in Section [1] the composite nature of the kinemat- 
ics in these multi-component systems can significantly affect 
the inferred rotational properties, and it is only by separat- 
ing the distinct kinematic components as we have done here 
that we can expect to get a reliable measure of the disc ro- 
tation speed and hence the galaxy's circular speed. We are 
therefore in a good position to place SO galaxies reliably on 
the Tully-Fisher relation. 

The other element we need for these SO galaxies are 
measurements of their absolute magnitudes. We carry out 
this analysis using the K-hand photometry introduced in 
Section [2] to minimise the impact of dust obscuration. There 
is a risk that the relatively shallow 2MASS exposures might 
miss some of the flux from the outer parts of the galaxy, 
so we use the magnitudes derived from the full GALFIT 
models obtained in Section [2] rather than the values quoted 
in the published 2MASS catalogue; in fact, this effect turns 
out to be rather small, with an average offset of only 0.05 
magnitudes. For NGC 3489, we had to use an SDSS 2;-band 
image to carry out the GALFIT modeling, so for this galaxy 
we calculate a total z-band magnitude from the model, then 
use the SDSS/2MASS integrated colours to convert it to a 
i^-band magnitude, which implicitly assumes that there are 
no strong colour gradients in this galaxy. 

To convert these values to absolute magnitudes, we 
use the distance modu li deriv ed from surface brightness 
fluctuations from Ton rv et all (|2001), systematically de- 
creased by 0.16 ma^ to take into a c count the updated 
C epheid zero poi nt of ^Free dman et all (12001), as discussed 
in ICoccato et al\ ([2009) and ICortesi et al\ (|2013h . Here, it 
is important to use a distance scale that has not itself been 
based on Tully-Fisher or other kinematic scaling relations, 
to avoid introducing a circularity into the analysis. 

Figure [6] shows the resulting Tully-Fisher relation for 
the SO galaxies analyzed here, compared to the relation fo r 
spiral galaxies (JTullv fc Piercell2000l : [kothberg et a/.ll2000l ). 
calibrated to the same absolute magnitude scale. As found 
previously by Ar agon- Salamanca (2008), there is clearly an 
offset between the SO and spiral relations, in the sense that 
the SOs are systematically fainter. This offset fits straightfor- 
wardly with a picture in which the SO galaxies ceased form- 
ing stars at some point in the past, leading to a steady fad- 
ing of their stellar populations. Again, there is no indication 
of any variation in this offset with environment, suggesting 
that this is not a driving factor in the timing of the tran- 
sition, but there does seem to be a trend in the sense that 
the most massive galaxies display the largest offset; in the 
context of the time since evolving off the spiral Tully-Fisher 
relation, this trend can be interpreted as another instance 
of "downsizing," in which more massive galaxies go through 
evolutionary transitions at an earlier stage than less massive 
galaxies. 
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Figure 7. X-band Faber-Jackson relation plotting spheroid ab- 
solute magnitude as a function of its dispersion for the SO 
galaxies studied here. Comparab l e kin ematic data for ellipti- 
cal galaxies from iMobasher et al\ (|l999l ) (filled red circles) and 
iMatkovic fc Guzman! (|2005l ) (open red circles) are also shown. 
The line shows a fit to the brighter elliptical galaxy data. 



^.S.2 Faher- Jackson relation 

Having separated spheroid and disc components, we are 
also in a position to look at connections in the other di- 
rection along the Hubble sequence to see how closely SOs 
might be related to elliptical galaxies. For these systems, 
the equivalent scaling is the Faber-Jackson relation between 
the spheroidal component velocity dispersion and its abso- 
lute magnitude. 

Accordingly we have plotted these quantities for the 
SO galaxy spheroids in Figure [71 using the spheroid disper- 
sions calculated in Section [4]T] and spheroid absolute magni- 
tudes derived as described above. For comparison, we have 
also plotte d the relation fo ll owed by elliptical galaxies as 
derived by iMobasher et all (|l999l ). These latter data are 
based on velocity dispersions measured at an effective ra- 
dius, away from any central spike in velocity dispersion, so 
should be directly comparable to the dispersions determined 
for the SO spheroids. These elliptical data do not span the 
full range of magnitudes seen for the SO spheroids, so we 
supplement them with kinematic data on fainter early-type 
galaxies in the Coma Cluster ( Matkovic fc Gu z man 2005) 
using K-band magnitudes from ISkrutskie et al\ ((2006). For 
these fainter galaxies, the kinematics were derived from fibre 
spectra, but the large three-arcsecond fibres used will cover 
most of these small faint galaxies at the distance of Coma, so 
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Figure 8. Plot showing how the SO galaxies are offset in mag- 
nitude from the mean Tully-Fisher and Faber-Jackson relations 
for spiral galaxies and elliptical galaxies respectively. 



should again be comparable to the luminosity- weighted aver- 
age values of asph derived for the SO spheroids. To obtain the 
absolute magnitude for the Coma galaxies, we placed the ap- 
parent magnitu d es tab ulated in Mobasher et al (1999) and 
ISkrutskie et all (|2006l l at a distance modulus of ??i — M = 
35.06, as also found using surface brightness fluctuations 
([Thomsen et a/.l [19971 1. The fact that the two comparison 
data sets follow a common relation even though their kine- 
matics were derived with very different radial weightings 
indicates that the results are unlikely to be sensitive to the 
exact manner in which the average velocity dispersion is ob- 
tained. 

The spheroids of the SOs follow a trend similar to 
the Faber-Jackson relation for ellipticals, although they 
seem to lie along its upper envelope. One possible expla- 
nation for such an offset might be that the influence of 
the disc surrounding the SO spheroid might in some way 
decrease its velocity dispersion. Howeve r, simulations by 
IPebattista. Kazantzidis fc van den BoschI (^2013. ) show that 
the formation of a disc around such a spheroidal component 
would in fact serve to compress it and increase its velocity 
dispersion, so the presence of a disc does not explain the dif- 
ference. A simpler explanation is that the offset is along the 
other axis, due to some enhancement of the spheroid light 
in SOs when compared to other elliptical systems. 



5 DISCUSSION 

In this paper, we have analysed the kinematics of a sam- 
ple of six lenticular galaxi es us i ng th e planetary nebulae 
data presented in Cortesi e t al\ ((2013|) to determine their 
dynamical properties out to large radii. This analysis has 
emphasised the importance of treating the spheroidal and 
disc components of these systems as distinct entities both 
photometrically and kinematically, and doing a careful job 
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of separating them. Through this process, the rather hetero- 
geneous properties of SOs over-ah start to be resolved into 
more consistent sub-components, and a reasonably coherent 
picture begins to emerge. 

The common factors we have uncovered are that: 

(i) The discs of SO galaxies are comparable to those of 
spirals, with similar flat rotation curves and falling veloc- 
ity dispersion profiles, but with a larger amount of random 
motions. 

(ii) The spheroids of SO galaxies show flat dispersion pro- 
files, similar to what is found in some ellipticals. 

(iii) SO galaxies follow the Tully-Fisher relation, but off- 
set to fainter magnitudes than spiral galaxies, with a greater 
offset for more massive galaxies. 

(iv) SO galaxy spheroids follow the Faber-Jackson rela- 
tion, somewhat offset to brighter magnitudes than elliptical 
galaxies. 

(v) There is no strong evidence that any of these effects 
depend on the current environment of the galaxy. 

The Tully-Fisher and Faber-Jackson findings are sum- 
marised in Figure [51 which shows the offset from each re- 
lation for the SO sample. In the absence of any systematic 
effect, we would expect these points to appear equally in all 
four quadrants, but in fact they appear strongly clustered 
in the region of the plot where the spheroids are too bright 
but the overall galaxies are too faint. 

With this level of detail starting to become available, 
we can begin to sketch out a plausible evolutionary sequence 
leading to the formation of SO galaxies. If these systems be- 
gan their lives as spiral galaxies, at some point they under- 
went a transition that removed their gas supply and hence 
cut off star formation. The process responsible does not de- 
pend strongly on current environment, but does seem to 
have acted earlier on more massive galaxies as another exam- 
ple of downsizing. The hotter discs of SOs could indicate that 
their high-redshift spiral progenitors had less ordered discs 
([Kassin et a/.ll2012l l , and that this larger amount of random 
motion became locked in when star formation ceased. How- 
ever, if the transformation occurred at a point when these 
systems had settled down to today's cold spiral systems, the 
hotter discs would then indicate that the galaxies had un- 
dergone a process sufficiently disruptive to heat the system 
significantly, but less violent than a merger. 

This aggressive process might well also dump some frac- 
tion of the gas being stripped from the disc into the centre 
of the galaxy, causing a last burst of turbulent star forma- 
tion there, enhancing the spheroid's luminosity significantly, 
shifting it off the Faber-Jackson relation. Such a scenario 
would also fit with the recent finding that the spheroids of 
SO galaxies have systematically young er stellar populations 
than their discs (Johnston et al. 201z). 

Of the various mechanisms advanced for causing galaxy 
transformation, perhaps the most plausible to realize such 
a scenario i s a f orm of the "harassment" invoked by 
iMoore et al\ (|l996l l as a process for transforming galaxies 
in clusters. In their paper, they demonstrated that this pro- 
cess of repeated high-speed encounters could strip gas from 
a galaxy while dumping some at the centre of the system 
under transformation as we require here. However, we also 
know that the process needs to be somewhat gentler than 
they modelled because we need to preserve the disc structure 



of the galaxy with only a fairly modest amount of heating. It 
remains to be seen whether such mild harassment, perhaps 
better thought of as "pestering," could be made to work 
across the range of environments in which SOs are found, 
but we do at least now have a growing number of obser- 
vational constraints against which any such model can be 
tested. 
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